
# coding: utf-8

# In[2]:


# -*- coding: utf-8 -*-
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D


# In[3]:


L = lambda F, D: 32.45 + 20*np.log10(F) + 20*np.log10(D)


# In[4]:


f = np.arange(3, 30)
d = np.arange(0.1, 1000)
f, d = np.meshgrid(f, d)
l = L(f, d)


# In[5]:


fig = plt.figure()
ax = Axes3D(fig)
ax.plot_surface(f, d, l, rstride=1, cstride=1, cmap='rainbow')
plt.show()

